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Abstract 

It is shown that the fluctuation in the temperature of the cosmic microwave 
background in any direction may be evaluated as an integral involving scalar 
and dipole form factors, which incorporate all relevant information about 
acoustic oscillations before the time of last scattering. A companion paper 
gives asymptotic expressions for the multipole coefficient Ci in terms of these 
form factors. Explicit expressions are given here for the form factors in a 
simplified hydrodynamic model for the evolution of perturbations. 
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I. INTRODUCTION 



The purpose of this paper is, first, to exhibit a general formahsm, ex- 
pressing the observed fluctuations in the cosmic microwave background tem- 
perature in terms of a pair of form factors, and then to carry out an illustra- 
tive approximate analytic calculation of these form factors. A companion 
paper [1] gives general asymptotic formulas for the coefficient Ci of the ^th 
multipole term in the temperature correlation function for arbitrary form 
factors, and also uses these formulas to calculate for the form factors 
found in the present paper. 

In Section II we show that under very general assumptions the fractional 
variation from the mean of the cosmic microwave background temperature 
observed in a direction n takes the form 

F{k) +in-kG{k)\ . (1) 

Here (Ia is the angular diameter distance of the surface of last scattering J**| 
and k^ek is proportional (with a k-independent proportionality coefficient) 
to the Fourier transform of the fractional perturbation in the total energy 
density at early times. (There are additional terms in AT/T that arise from 
times near the present, and chiefly effect the multipole coefficients Ci for 
small i, especially ^ = and £ = 1. These will be discussed in Section IV 
and in the Appendix. Effects from a changing gravitational field soon after 
the time of last scattering are included in Eq. (1).) 

One advantage of this formalism is that it provides a nice separation 
between the three different kinds of effect that influence the observed tem- 
perature fluctuation, that arise in three different eras: (1) at very early 
times, (2) during the era of acoustic oscillations, and (3) from the time of 
last scattering to the present: 

1. The k-dependence of the unprocessed fluctuation amplitude ek reflects 
the space-dependence of fluctuations in the energy density at very early 
times. The average of the product of two es is assumed to satisfy the 
conditions of statistical homogeneity and isotropy: 

^kCk') = 5'(k + k')P(fc) (2) 

** Note that in speaking of a surface of last scattering, we are not necessarily assuming 
that the transition from opacity to transparency takes place instantaneously. The physical 
wave number vector k varies with time as l/ffl(f) (where a{t) is the Robertson-Walker scale 
factor), while for large redshifts dA varies as a{t), so the product cLa^ is nearly independent 
of what we choose as a nominal time of last scattering. 
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with k = |k|. Since the reahty of the fluctuations in the energy density 
requires that ej^ = e_k, the power spectral function V{k) is real and 
positive. It is common to assume a "straight" spectrum 

V{k) oc yt"-^ (3) 

For instance, the 'scale-invariant' n = 1 form[2] suggested by theories 
of new inflation[3] is: 

V{k)=Bk-^, (4) 

with B a constant that must be taken from observations of the cosmic 
microwave background or condensed object mass distributions, or from 
detailed theories of inflation. 

2. The form factors F{k) and G{k) characterize acoustic oscillations, with 
F{k) arising from the Sachs-Wolfe effect and intrinsic temperature 
fluctuations, and G{k) arising from the Doppler effect. 

3. Taking account of a constant vacuum energy density as well as cold 
matter in the time after last scattering, it is easy to calculate the 
angular diameter distance of the surface of last scattering: 



dA = —rj-r, — sinh 



dx 



(5) 

where ~ 1100 is the redshift of last scattering, = 1 — — ^Mj 
and r^A and Q^m are as usual the present ratios of the energy densities 
of the vacuum and matter to the critical density 3i?o / SttG. 

We note in particular that F{)z) and G{k) depend on and on the 

baryon density parameter ri^/i^ (where h is the Hubble constant in units of 
100 km/sec/Mpc), but since the curvature and vacuum energy were negligi- 
ble at and before last scattering, Fik) and Gik) are essentially independent 
of the present curvature and of JIa- The exponent n in is expected to 
be independent of all these parameters. On the other hand, dA is affected by 
whatever governed the paths of light rays since the time of last scattering, 
so it depends on ^^A, and the curvature parameter ilc", but it is essen- 
tially independent of quantities like the baryon density parameter that 
effect acoustic oscillations before the time of last scattering. In quintessence 
theories dA would be given by a different formula, but 7^(A;) and the form 
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factors would be essentially unchanged as long as the quintessence energy 
density is a small part of the total energy density at and before the time of 
last scattering. 

Another advantage of this formalism is that, although Ce must be cal- 
culated by a numerical integration, it is possible to give approximate ana- 
lytic expressions for the form factors in terms of elementary functions. The 
detailed confrontation of observation and theory must necessarily be done 
using computer codes that take into account all relevant astrophysical and 
observational effects. [4] Nevertheless, there is some value in also having an 
analytic treatment that, though not as accurate as possible, is as simple as 
possible while still capturing the main features of what is going on. The 
point is not to compete with the computer codes, but rather to gain some 
feeling for what is going on, in order to help us judge how predictions for 
the cosmic microwave background fluctuations may change with alterations 
in the underlying assumptions. 

Analytic treatments of fluctuations in the cosmic microwave background 
already exist in the literature [5]. Our main purpose in going over the same 
ground here is not to give a more accurate or comprehensive treatment of 
acoustic oscillations, but to obtain simple expressions for the form factors 
as examples to which to apply the asymptotic formulas for derived in 
reference [1]. To derive analytic expressions for the temperature fluctuation 
it is necessary to neglect the contribution of radiation and neutrinos to 
the gravitational field, which should be a fair approximation near the first 
Doppler peak but not much beyond that. We employ a purely hydrodynamic 
treatment, relying on the Boltzmann equation only implicitly in the values 
used for the shear viscosity and heat conduction coefficients; the effects of 
viscosity and heat conduction are included from the beginning, not just by 
inserting damping factors; and 'Landau' damping due to the finite duration 
of the era of last scattering is included along with 'Silk' damping due to 
shear viscosity and heat conduction. As far as I know, this is the first work 
to obtain explicit analytic expressions for the temperature fluctuations that 
are correct within these approximations. 

Section III presents an analytic calculation of the evolution of pertur- 
bations in synchronous gauge up to the time of last scattering, which is 
then used in Section IV to calculate the form factors. For very small wave 
numbers the form factors are found to be 

F(k) ^ l-3kHl/2-3[-r' + r^Hl + 0]k%/i + ... , (6) 
G{k) ^ 3ktL - 3kHl/2{l + e) + . . . , (7) 
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while for wave numbers large enough to allow the use of the WKB approxi- 
mation, i. e., Hl > ^, the form factors are[6] 



-3C + + (1 + 0-'/^e-'=''^i^ cosikdn) , 

(8) 
and 

G{k) = (1 + O'^^H'^ + 2i/kHl)-^e-^'''^o sm{kdH) ■ (9) 

Here tL is the time of last scattering; ^ is 3/4 the ratio of the baryon to 
photon energy densities at this time: 

C=f^') ^27nnh'; (10) 

dn is the acoustic horizon size at this time, given by Eq. (75), and dc is a 
damping length, given by Eq. (89). 

II. FORM FACTORS 

We first justify Eq. (1) under very general assumptions, not limited to 
those of Section III. At and before the time of last scattering the spatial 
curvature was negligible, so small perturbations in the cosmic metric and 
in all particle distributions at these times may conveniently be expressed 
as Fourier transforms of functions of a co-moving wave number vector q 
and the time t. Effects like pressure forces that involve spatial gradients 
are important for a given q only when the physical wave number q/a{t) is 
at least as large as the cosmic expansion rate, which is of order 1/t. Since 
a{t) vanishes for t — no more rapidly than -y/i, the ratio qt/a{t) vanishes 
as t — > 0, so whatever the value of q, there will always be some time early 
enough so that pressure forces and other effects of spatial gradients are neg- 
ligble. At such early times, perturbations grow or decay with powers of time. 
Generically there is one most rapidly growing mode, and this is the one that 
eventually grows into the perturbations seen at the time of last scattering. 
Since the equations for the time dependence of the perturbations are linear, 
the Fourier transforms of all perturbations to the metric and particle dis- 
tributions during the era of last scattering will then be proportional to the 
Fourier transform Cq of any one of these perturbations at any sufficiently 
early time. For definiteness, we can take Cq to be q^^ times the Fourier 
transform of the fractional perturbation to the total energy density at some 
very early time, a choice that will prove' to be convenient in Section IV. 



F{k) = il + 2^/kHl)-' 
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Since the fractional change in the observed microwave background tem- 
perature seen in a direction n is linear in the perturbations to the metric 
and photon and matter distributions at various times during the era of last 
scattering, it can be written as 



where r(t) is the co-moving radial coordinate of a source scattering hght at 
time t that would be received at the present. Note that the quantity J can 
depend on q only through the scalars q and q ■ h, because the differential 
equations governing the growth of perturbations are rotationally invariant, 
even though the initial fluctuation amplitude Cq is not. 

We can make a great simplification in Eq. (11) by taking advantage of 
the fact that the Robertson-Walker radial coordinate r{t) is nearly constant 
during the era of last scattering. Using equilibrium statistical mechanics to 
calculate the hydrogen ionization, and simple Thomson scattering to calcu- 
late scattering probabilities, one finds that for Q,b/^m = 0.2, the probability 
that a photon will never again be scattered rises from 2% at 3360°K to 98% 
at 2780°K. (This depends on the assumed value of JIb/Om, but very weakly; 
for instance, for Qb/^m = 0.12, the probability of no future scattering rises 
from 2% to 98% as the temperature drops from 3400°K to 2810°K.) For 
definiteness, we will round off these temperatures, taking the era of last 
scattering to extend from a temperature of 3400° K down to 2800° K, corre- 
sponding to a drop in rcdshift z from 1220 to 1010. The radial coordinate 
can be expressed in terms of z by the well-known formula 



where Clc = 1 — Oa — JIm, a{t) is the Robertson-Walker scale factor, tz is 
the time corresponding to redshift z, and to is the present. This approaches 
a constant limit for z ^ oo, and therefore varies very little in the range 
from z = 1010 to z = 1220 (or, for that matter, even in the range from 
z = 1010 to z ^ oo.) For instance, if we take the popular values Q,m = 0.3 
and fly = 0.7, then the fractional change in r{tz) as z drops from 1220 to 
1010 is 0.0034. We can therefore take the exponential in Eq. (11) outside 
the time integral, so that 




(11) 




(12) 




(13) 
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where is any conveniently chosen time during the era of last scattering, 
say at a redshift zl = 1100. 

It is convenient to replace the co-moving wave number vector q with the 
physical wave number at last scattering: 

k = q/a(ti) . (14) 

Eq. (13) may then be written 

= J d^k e^^-^'^^ ek J dt J{k, k-n,t), (15) 

where J^{k, k -h^t) = a{tL)^J{q, q ■ h, t), = Cq, and dA = r{ti,)a{tL) is the 
angular diameter distance of the surface of last scattering. 

Part of the observed temperature fluctuations arise from perturbations 
in scalar quantities, like the gravitational potential and the intrinsic tem- 
perature, and therefore makes a contribution to J that is independent of 
n. Another part arises from fluctuations in a vector, the velocity of the 
baryon-electron plasma, and therefore makes a contribution that is linear 
in n. Leaving aside other effects like gravitational radiation, the function J 
therefore takes the form 

J{k, k-n,t)= r{k, t)+ik-n g{k, t) , (16) 

with ^{k, t) arising from the Sachs- Wolfe effect and intrinsic temperature 
fluctuations, and Q{k,t) arising from the Doppler effect. Using this in 
Eq. (15) then gives 

= J d^k e*-^*^^ [F{k) +ik-h G{k)\ ek , (17) 

which is the same as Eq. (1), with the form factors identified as time integrals 

F{k) = J dt J^{k,t) , G{k) = J dt g{k,t) . (18) 

This time integration introduces a damping of oscillatory part of the form 
factors, but this will be less important than the effects of heat conduction 
and viscosity in the time interval between recombination and last scattering. 

III. EVOLUTION OF PERTURBATIONS IN SYNCHRONOUS 

GAUGE 
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We now turn to the approximate analytic calculation of the form factors. 

A. General Approximations 

We make two assumptions that will allow great simplifications in this 
calculation: 

1. The contents of the universe up to the time of last scattering are taken 
to consist of collisionless cold dark matter, collisionlcss neutrinos, a 
baryon-electron plasma treated as a perfect fluid, and a photon gas 
coupled to the plasma by Thomson scattering, with a short but non- 
negligible photon mean free time. The finite duration of the era of 
last scattering, when the mean free time becomes too large to allow 
a hydrodynamic treatment, will be taken into account by the time 
integrals in Eq. (18). 

2. It is assumed that only the cold dark matter contributes to the ex- 
pansion rate of the universe before the time of last scattering and to 
perturbations in the metric. This is not a very good approximation, 
but it is the price that has to be paid to get analytic expressions for 
the observed temperature fluctuation. To minimize errors introduced 
by the incorrect treatment of acoustic oscillations before the cross-over 
time tc when the photon energy density equaled the dark matter en- 
ergy density, it is necessary to restrict the wave number to be less than 
an upper bound given in Section V. 

B. Gravitational Field 

We begin by reminding the reader of the equations that govern pertur- 
bations in the metric and fluid properties before the time of last scattering. 
The perturbed metric is taken as 

Qixv total (x, t) = gfj,t,{t) + hf^i^ipc, t) , (19) 

where g^^ is the Robertson-Walker metric in co-moving coordinates x with 
spatial curvature neglected: 

500 = -1 , ffoi = , g,j = a^it)5ij , (20) 

and /i^jy(x, t) is a small perturbation. We work in synchronous gauge, defined 
by the conditions 

hoi = hoo = Q, (21) 
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and by the requirement that the cold dark matter particles have time- 
independent spatial coordinates. These conditions leave an unbroken resid- 
ual gauge invariance, under the transformation 

with ej an arbitrary function of x but independent of t. As we will see, in 
synchronous gauge the evolution of the compressional modes that concern 
us here depends on the gravitational field only through a quantity that is 
invariant under these these transformations 



, d fhkk 



(23) 



The spatial curvature is negligible at and before the time of last scattering, 
so it will be convenient to express i/'(x, t) as a Fourier transform: 

V^(x,i) = / d^e'^-^i^^it). (24) 

Likewise, the total proper energy density of each of the constituents of the 
universe (labelled f = D, B, ^ for dark matter, the baryon-electron plasma, 
and photons, respectively) is written 

Qf totai(x, t) = e/(t)+%(x, t) , (5p/(x, t) = j (fq gfq{t) e^^"' , (25) 

with quantities carrying a subscript q denoting small perturbations. Under 
our second assumption, the gravitational field equation in synchronous gauge 
reads [7] 

^ (a^V'q) = -47rGaVDq • (26) 
C. Dark Matter Perturbations 

The dark matter particles are assumed to ride on the expanding co- 
ordinate mesh, with negligble peculiar velocities. (This is not affected by 
perturbations to the gravitational field, because in synchronous gauge these 
perturbations leave Fqq zero.) Hence their energy- momentum tensor has 
only a 00-component, = qd totsd- For the metric (19)-(21), the energy 
conservation equation TjJ\^ = then reads 
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or in other words 

-V'q , (27) 



where 6dci fractional dark matter density perturbation: 

feq = PDq/PD ■ (28) 

Combining Eqs. (26) and (27) gives 

(a2^)=47rGaWDq. (29) 



dt 



During the dark matter dominated era, a oc i^/^ and inGpo = 2/3i^, so 
Eq. (29) can be written 

As is well known, the two solutions go as and i^/^. If these two modes 
have comparable strengths for very small t, then the relevant solution is the 
one that is most rapidly growing, which we shall write as 

Sdci = ^q*'^' ' V'q = - ^iVqi- V3 . (31) 

(The normalization constant Nq will play a role in this section similar to 
that of the constant Cq in Section II.) 

D. Plasma and Photon Perturbations 

Next, let us consider the imperfect fluid formed by the baryon- electron 
plasma and the photons. It has a total velocity four- vector of the form 

KtJ^^ t) = U^ + J d\ U^{t) e^^-- (32) 

where is the unperturbed velocity four-vector 

[/° = 1 , W = 0, (33) 

and Uq^{t) is a small perturbation. The normalization condition g^j, total f^totai^totai 
— 1 tells us that the first-order perturbations are purely spatial 

C/°(t) = . (34) 



9 



We will be considering only compressional modes, so we will assume that 
the spatial part of U^otai gradient of a velocity potential u: 



(35) 



We will write the conservation laws for this fluid in terms of fractional pertur- 
bations to the baryon-electron plasma mass density and the photon energy 
density: 

= PBq/PB , Sjq = Pjq/ Pj ■ (36) 

The particle conservation equation [8] for the baryon-electron plasma mass 
density is then 



dt 



q Uq-tpd- 



(37) 



The energy conservation equation [9] for the baryon-electron-photon fluid is 



d_ 
di 



(^PbSbci + Pi^iq) + ^ (^PsSBq + ^P7'^7q^ = " (^PB + 



->pl 



T d{a?Uc 
dt 



+ 



T5^ 
4a2 



(38) 



where T is the unperturbed photon temperature and x is the coefficient of 
heat conduction caused by photon energy transport. Finally, the momentum 
conservation equation [10] is 



4 + 167rGr/ 
dt 



-a [pB + '^Pj 



d_ 
di 



_ 1 3 

— P^'^^q 



-q Mq + Ipq 



(39) 



where rj is the coefficient of viscosity due to photon momentum transport. 
By using Eq. (37) and recalling that pB oc a~^, oc a~^, and T oc a~^, we 
can simplify Eq. (38) to read 



d 






di 




Pi 



■^|(«(^Bq+^)) 



4a2 



(40) 



Also, using Eq. (37) lets us write Eq. (39) as 



d_ 
'di 



+ levrGry 



a \PB + -Pi 



d_ 

dt' 



q^a? 



P-ySjq + ^V^Bq ■ (41) 
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Now, rj/p^ and / p-y are of the order of the photon mean free time, 
which as long as hydrodynamics is appUcable must be short compared with 
the cosmic age. Therefore we can neglect 77 and x everywhere, except where 
they are accompanied with a maximum number of space and /or time deriva- 
tives of 6bci '-'^ *^7q' which case powers of a high wave number can com- 
pensate for the smallness of x or V- Then Eqs. (40) and (41) simplify further 
to 

d_ 
di 



-Sbci ~ 



(42) 



d_ 
dt 



(pB + ^Pij (5sq + Ipd) 



-pry 5, 



7q 



xTa' 




7 , 



(43) 



We also neglect terms of second order in x and/or r/, so we can set 
equal to 3(^-yq/4 in the dissipative terms in Eqs. (42) and (43). Then using 
Eq. (42) to eliminate (5sq in Eq. (43) gives our differential equation for 5^q: 



d_ 
'di 



PB + -P7 



)( 



3d(5. 



7q 



4 dt 



V'q) 

d5. 



4P7 



4 dt^ 4a2 dt 



2„3 "'"7q 

dt 



(44) 



It will be convenient to multiply with t^^^/a^p^. Recalling that p^ oc a ^, 
PB oc and a oc i^/^, this gives finally 



dt 



(l+i.).2/3^ 



7q 



+ 



Pj 



dt 



+ 

4*2/3 d 



f 0.4/3 d^^-yq , ,2.4/3 '^'^7q ) 
4/9^ dt^ dt J 



8N, 



^{1 + 3R) 



(45) 



3 dt\' ' 27 
where (as before) ti is some typical time of last scattering, k = q/a{tL) = 
qt^l^/t^l^a{t), and 

i? = 3/9ij/4p^ oc a . (46) 

We next turn to two different ranges of wave number in which it is 
possible to find an analytic solution of this equation. 
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E. Solution for large k 



We consider first wave numbers that are large enough to allow the use 
of the WKB approximation. For this purpose, we introduce a new variable 

C^{t/tL)'/'. (47) 

(This is the usual conformal time rj, but with a different normalization.) 
Multiplying Eq. (45) with 9t^^ then gives 



A. 



'4^. W dc^ dc 



P-rC dC 

+ 3fcV^)=^(l + 3eC^). (48) 



Here we have again assumed that dissipative terms are negligible except 
where a maximum number of derivatives (i.e., factors of k and/or (^-derivatives) 
acts on 5yq. We have also used the fact that R oc a to set R = , where ^ 
is the ratio (46) at time tL. 

In the absence of dissipation, Eq. (48) would have the exact solution 

^ _ 8iVqtf(l + 3ga 

9(Pti + 20 ■ ^^^^ 

(This is actually independent of our choice of t^, because ^ and k'^t'j^ both 
scale as .) The neglect of dissipation is justified in this solution, because 
the rate of change of this expression does not yield a factor of the large wave 
number k that could compensate for the smallness of x and rj. 

To this particular solution, we must add a suitable solution of the cor- 
responding homogeneous equation. In the absence of damping we can find 
exact solutions of the form Pi/(i\/fC)) where Pu{z) is the usual Legendre 
function, and u is either of the roots of the quadratic equation i/(i/ -|- 1) = 
— 'ik'^t\/^. But this will not be useful in calculating the in our companion 
paper [1]. To get a more useful result, we must use the WKB approximation. 

Under the assumption that 

ktL > i (50) 
we can find a pair of approximate solutions of the homogeneous equation 

(5-yq oc exp(±i(/3) (51) 
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with 



Jo 



provided we neglect dissipative terms. Note that if as wch as r] and x were 
zero, then these homogeneous solutions would be exact. More generally, 
inspection of Eq. (48) shows that these are approximate solutions if the 
fractional rate of change of 1 is small compared with the rate of change 
of the phase (p: 

2gC < ktLV3 
which is true at all times if and only if it is satisfied at ( = 1, i.e.. 



ktL > 



2^ 



For plausible values of ^ this condition is actually somewhat weaker than 
Eq. (50), but we will need the greater strength of Eq. (50) later, when we 
calculate the plasma velocity potential. 

We can do better than Eq. (51), and include the effects of viscosity and 
heat conduction, by seeking solutions of the form 6-yq = Aexp{±iip), with 
A a slowly varying real amplitude. By calculating the rate of change of the 
Wronskian of these two solutions (and replacing d^(5^q/(iC^ with —3k^t'j^{l + 
^C'^)~^dS^ci/dC in the dissipative term), we easily find the WKB solutions of 
the homogenous equation: 



^^qOc(l+^C')"'/'exp 



where 



= 3tL 



+ 



2^2 



The viscosity and heat conduction coefficients are given by [11] 



16 
45 



(53) 



(54) 



(55) 



Here r-y is the photon mean free time 



{27rmekBTo/h''r/\l + ZLr/\l-Y/2y/^ ^ [2kBTL J ^ 
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where gt is the Thomson scattering cross section, h is (only here) the orig- 
inal Planck constant, Tq = 2.738°K is the present microwave background 
temperature, Tq = 3100°K is the temperature at last scattering, ks is Boltz- 
mann's constant, Y ~ .23 is the primordial helium abundance, and A = 13.6 
eV is the hydrogen ionization energy. 

The relevant solution is again the one that increases most rapidly at 
early times, which we can find by requiring that — > as ^ — ^ 0. In the 
limit C ^ the phase (p vanishes as 0(C), while vanishes more rapidly 
because the mean free time of photons is very small at early times. Hence 
the linear combination of the particular inhomogeneous solution (49) and 
the homogeneous solutions (53) that grows most rapidly at early times is 



9(A;2ti + 20 



l + 3eC^-(l + eC')-VV=^^%os^ 



(57) 



(We would be able to neglect the term 2^ in the denominator only under 
the condition kti > -\/2^) which for plausible values of ^ is stronger than 
our assumption (50). ) 

To calculate the velocity potential of the plasma-photon fluid for large 
wave numbers, we will also need the rate of change of c^sq. At times of order 
tij the time derivatives of ^C^i V') ^re of the orders of ^/II, k, and 

tl, respectively, where tl is the photon mean free time t ^ rj/pj ^ xT/Pj 
at time ti. We arc assuming that ktL > so the time derivative of ip is 
larger than the time derivative of^(^. Eq. (56) shows that damping becomes 
important if k'^tiTi > 1, , but even for such large values of k we can still 
limit ourselves to the case 

fcTL < 1 , (58) 

in which case the time derivative of ip is also larger than the time derivative 
of k'^V'^. Hence for wave numbers k in the range defined by Eqs. (50) and 
(58), we have 

d5^ ^ 8iVqtf fce-^^^%iny 

dt Q^{i+ief/\kHi+2oe' ^ ' 

The dissipative terms in Eq. (42) are smaller than this by a factor fcr, so 
here wc can take ^gq — 3(j-yq/4, and Eqs. (37), (31) and (59) then give the 
velocity potential 



2Arq 



3Pa2(t^)t^/\ 



-1 + ^ 



V3(l + eC')3/4(A:2t2 +20C. 



(60) 
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F. Solution for small k 



Here we can neglect viscosity and heat conduction. For k = 0, Eq. (48) 
has an obvious solution 



= 4iVqi^/\V3 = 4feq/3 . 



(61) 



To this we can add any linear combination of the two solutions of the 
corresponding homogeneous solution, for which Sjq is respectively time- 
independent or proportional to 



1 + ' 

which near the beginning of the dark-matter dominated era goes as (. As 
^ increases these homogeneous solutions become negligible compared with 
the inhomogeneous solution (61), so at later times the solution for = is 
given by Eq. (61). 

To get the term in Sjq of first order in k^, we can use the solution (61) 
in the terms in Eq. (48) proportional to k^, so that 



d_ 



skHlSod 



(62) 



Discarding a homogeneous term for the same reason as before, we have 



4-j-s. 



which gives 



dC V4 



kHl 



1 + 



in(i + ec') +... 



(63) 



(64) 



Also, Eqs. (37), (27), (63) and (31) give the plasma velocity potential for 
k ^ as 

' (Sb^ - Sud) - -^37^^ • (65) 



^ q^dt ""'^^ • 3a2(ti)(l+^C') 

As we will see, this provides a small correction to the Doppler shift, which for 
small A; will turn out to be mostly due to perturbations in the gravitational 
field. 
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IV. OBSERVED TEMPERATURE FLUCTUATIONS 



There are three separate sources of the observed temperature fluctuation 
in the cosmic microwave background: the Sachs- Wolfe effect due to pertur- 
bations in the gravitational potential, the Doppler effect due to plasma pe- 
culiar velocities, and the intrinsic temperature fluctuations themselves. We 
will consider each of these in turn, and then put the results together. In 
calculating the Sachs-Wolfe and Doppler contributions, we will use a non- 
relativistic approach, taking the effect of the gravitational field perturbations 
on the observed photon temperature to consist entirely of the time dilation 
caused by a Newtonian gravitational potential plus the Doppler shift caused 
by the gravitational acceleration of the source and receiver. This approach 
has the virtue of getting useful results quickly, but the results obtained in 
this need to be justified by a thoroughly relativistic treatment of the Sachs- 
Wolfe and Doppler effects, which will be given in the Appendix. 



We can define a Newtonian gravitational potential (p as the solution of 
the Poisson equation 



with the factor inserted to take account of the difference between the 
Robertson- Walker co-moving coordinate vector x used here and the coordi- 
nate vector a(i)x that measures proper distances at time t. Using Eqs. (31) 
and (25), this gives 



It is important to note that this is time-independent during the dark matter 
era, when a{t) oc i^/^. 

This potential makes two separate contributions to the Sachs-Wolfe ef- 
fect. There is a gravitational redshift, yielding a fractional fluctuation in 

the observed temperature in a direction n equal to (j){rLh) — (/>(0), where 
is the Robertson- Walker radial coordinate of the surface of last scattering. 
There is also a time-delay; if the unperturbed cosmic temperature reaches 
the value Tl — 3000°K of last scattering at a time i^, then the gravitational 



A. Sachs-Wolfe Effect 



a-2(t)V2(/)(x,t) = AnG5pD{^,t) 



(66) 



0(x,i) 




(67) 
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potential causes the cosmic temperature in a direction h to reach the value 
Tl at a time [1 + (t){rLn) — (/)(0)]iL, so that the redshifted temperature seen 
now is changed by a fractional amount [12]: 



tLa{tL)/a{tL) <t>{rLh) - (/-(O) 



(This argument is valid only because is time-independent; otherwise we 
would have to consider the complete gravitationally delayed time-history 
of the cosmic temperature, as done in the Appendix.) Combining the two 
effects, the net fractional change in observed temperature is 



/Ar(n) 



1 



Sachs -Wolfe 



(t>{rLn) - <^(0) 



(68) 



As we shall see in the Appendix, this formula can be derived using the 
formalism of general relativity, which in synchronous gauge gives the famous 
factor of 1 /3 directly, without having to consider separately the gravitational 
red shift and time delay. 

It will be convenient to rewrite Eq. (67) in terms of the physical wave 
number at the time of last scattering, k = q/a{tL), so that Eq. (68) gives 



/AT(n) 



T 



Sachs-Wolfe 



J'k-hdA 



- 1 



(69) 



where ek is an amplitude for fluctuations not processed by acoustic oscilla- 
tions, defined by 



ek d'^k = 



and dA = rLa{ti) is the angular diameter distance of the surface of last 
scattering. 

B. Doppler Shifts 

The plasma velocity potential Uq calculated in Section III yields a pressure- 
induced plasma velocity perturbation 



' pressure 



(x,t) = a{t) V j d^q e''i '^Uq(t) 



(71) 



(The factor a{t) enters because it is the velocity in co-moving coordinates 
that is given by the co-moving gradient of the velocity potential.) This yields 
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a Doppler shift of the temperature of the cosmic microwave background seen 
in a direction n: 



(^^) = -ia{tL) fd^qh-ciu^{tL)e'^-'''^ 

V -1 / pressure Doppler J 

= -ijd^kfi-key, g{k) t^^'^'^^ , (72) 
with the form factor g(k) given by Eq. (65) for small k as 

gi^k) = (73) 

2(1 + ^) ^ ^ 

and for large k by Eq. (60) as 

g{k) = SktL - V3(l + + 2^/khl)-^e-''''^n sin (kdn) (74) 

where ^ as before is 3/4 the ratio of baryon and photon energy densities at 

the time of last scattering, is the damping length D given in Eq. (56), 
evaluated at (" = 1 (actually, as discussed in the next section, at ( a little less 
than unity) , and dn is the acoustic horizon at the time of last scattering: 

In the non-relativistic approach used here, there is also an additional ve- 
locity perturbation induced by the gravitational potential 0(x). The proper 
peculiar velocity Vgrav produced in this way is given by the equation of 
motion [13] 

|vg,,.(x,t) + ^vg,,.(x,t) = -^V<A(x) (76) 

Because the gravitational potential (p is time-independent, this has the sim- 
ple solution 

Vgrav(x,0 = -a-i(t)i V<^(x) = |^ J d\ q-^e'^'^N^ . (77) 

This contributes a fractional temperature shift seen in a direction n: 

/AT(n)~ 



V T 



= • Wgra-vinrL, II) - Vgrav(0, to)] 

gravity Doppler 

1/3, 



= 3i J d^k {k ■ h) ktL€]^ 



, (78) 
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where to is the present time. (A general relativistic derivation of this result 
in synchronous gauge is given in the Appendix.) 



C. Intrinsic Temperature Fluctuations 



The fractional change in the photon temperature is one-fourth the frac- 
tional change in the photon energy density. The contribution of intrinsic 
density fluctuations at the time ti of last scattering to the fractional change 
of temperature seen coming from a direction n is therefor^ 



^Ar(n) 



5p-y{nrL,tL) 



/ intrinsic 

Eqs. (57), (64) and (70) then give 
fAT{h 



- / (Pq e 



with the partial form factor / given by 



d^ke^fik)e 



S,^{tL) . (79) 



(80) 



fik) 



-3kHl/2 + 3 [r ' - ln(l + 0] + ... 

(1 + 2i/kHl)-^ f-1 - 3C + (1 + O-i/^e-'^'^i cos{kdH] 



D. Total Temperature Fluctuations 



k^Q 
k large . 
(81) 



We now put together the fractional temperature fluctuations given by Eqs. (69), 
(72), (78), and (80), and obtain the total fractional temperature fluctuation 



Ar(n) 
T 



d^k€k { F{k) +ik-hG{k) 



J'k-hdA 



1 - 3ik ■ 



n 



.1/3 



(82) 



t There is a subtlety here. To the extent that the opacity drops sharply from 100% 
to zero, last scattering occurs at a fixed value of the perturbed temperature T + ST, near 
3000°K, rather than at a fixed value of the unperturbed temperature or the time. The 
effect of the intrinsic temperature fluctuation 5T{t) is thus to change the time of last 
scattering, in such a way as to produce a change —ST in the value of the unperturbed 
temperature T{t) at this time. Since T{t) oc l/a{t), we then have Sa/a = +ST/T at the 
time of last scattering, so that the observed temperature is shifted by the change in the 
cosmological redshift by a fractional amount AT/T — Sa/a = +ST/T. 
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where F{k) is the total scalar form factor, given by Eqs. (69), (80), and (81) 
as 



l + f{k) 



1 



3[-r'+r'ln(l+0] kHi/A+... 



[l + 2i/kHi 



^ + 2i/kHl + {l + 0-^l^e 



cos{kdH] 



k^O 
k large . 

(83) 



and G{k) is the total dipole form factor, given by Eqs. (72), (73), (74), and 
(78) as 



G{k) = 3ktL-9ik) 



3ktL - 3k^tl/2{l + ^) + . . . k^O 
V3{1 + 0"^/^(l + 2i/kHl)-^e-^^T^l sm{kdH) k large. 

(84) 

The last two terms in the curly brackets in Eq. (82) contribute only to the 
multipole coefficients Q for £ = and I = 1[15], and may therefore be 
dropped (as they are in Eq. (1)) in considering the higher multipoles. 

We see that the WKB solution for large k gives a poor picture of what 
happens for A; — 0, except in the case ^ ^ 1, where dn = V^tL, in which 
case the above pairs of expressions for F{k) and G{k) agree for small k. 

As discussed in Section II, it still remains to average over the time of last 
scattering. The effect of this averaging on the damping factor exp — 
is small[15]. Otherwise, the averaging over t chiefly affects the sin kdn and 
cos kdn factors in Eqs. (83) and (84), which oscillate rapidly with the time 
of last scattering when k is large. We will approximate the probability 
distribution of the actual time of last scattering i as a Gaussian of the form 
(l/vr At) exp(— (t — tx)^/At^), where ti is a nominal time of last scattering. 
Replacing ti in the sines and cosines in Eqs. (83) and (84) with t, multiplying 
with this probability distribution, and integrating over t then gives the same 
result for the form factors for large k , but with an additional term now added 



to Vl: 



4 



At 
2t^ 



(85) 



This is a sort of 'Landau damping,' except that the damping arises from 
a spread in the time at which the temperature of the medium is observed 
rather than from a spread in wave numbers. As we will see in the next 
section, this term makes a smaller but not insignificant contribution to the 
total damping. 
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V. DISCUSSION 



In a companion paper[l] we show how to use the formula (82) for the 
total temperature fluctuation to derive expressions for the coefficient of 
the term of multipole number I in the temperature fluctuation correlation 
function for general form factors F[k) and G{k). As we will see there, the 
contribution of the scalar form factor F{k) to arises mostly from wave 
numbers of order i/dA (where cIa is the angular diameter distance of the 
surface of last scattering), while this approximation is much worse for the 
contribution of the dipole form factor G{k). 

For the present, we will content ourselves with noting that if we tenta- 
tively use the WKB approximation, neglect damping effects, and drop the 
terms in the second line of Eq. (83) proportional to ^//c^t|, then for ^ less 
than 0.311 (that is, for 3^ < (1 + the squared scalar form factor 

F^{k) has peaks at the wave numbers 

kn = n-K/dn , (86) 

(with n = 1,2, . . .), with higher peaks for odd n (where the two terms in 
F{k) have the same sign) than for even n. The minima are at the zeroes 
of F{k). For ^ > .311 the only peaks are those for n odd, and the minima 
are at n even. This suggests that there should be peaks in Q near = 
(2n — l)7rdA/dH and either lower peaks or dips near 2TnrdA/dH, depending 
on the value of ^. These peaks are known as the Doppler peaks (though 
Eq. (84) shows that the contribution of the Doppler shift is very small at 
all the wave numbers fc„.) These results depend critically on the negative 
sign of the term —3^ in the second line of Eq. (83); if this term had turned 
out to be positive then for ^ > .311 the positions of the peaks and dips 
would be interchanged. Despite what is sometimes said [16], there is no way 
without detailed calculations to see that the first Doppler peak should be at 
i 2± TrdA/dn, rather than at a multipole number twice as large. 

We can now check whether the WKB approximation used in subsection 
3E is valid at the first Doppler peak. According to Eqs. (86) and (75), we 
have 

kitL = , , (87) 

^/3in(v?+^/rTe) ^ ' 

so the ratio of the wave number at the first Doppler peak to the mininum 
wave number kmin allowed by the inequality (50) is 

" (88) 



A;min V3?ln(V! + v/rT^) ■ 
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The WKB approximation is valid at wave numbers down to the first Doppler 
peak if this ratio is sufficiently larger than unity. For instance, for Q,Bh^ = 
0.03 we have ^ = 0.81, so Eq. (88) gives ki/kmin = 2.5, making the WKB ap- 
proximation fairly good at the first Doppler peak. The WKB approximation 
is somewhat better at ki for smaller values of O^/i^, though it still breaks 
down at smaller wave numbers unless flBh? = 0. For all plausible values of 
^ the WKB approximation is excellent at the higher Doppler peaks. 

Next, let us consider the importance of damping. It might seem that we 
should calculate the damping length Vl by integrating in Eq. (54) up to the 
time of last scattering, corresponding to ( = 1. But at the nominal time 
of last scattering (defined so that the probability of any future scattering is 
50%), the photon collision rate l/r^ given by Eq. (56) is {).2yJQ.B/^M /tl^ 
which is already considerably smaller than the expansion rate 2/3tL, so that 
we cannot trust the hydrodynamic calculations used to obtain Eq. (54) . We 
will instead integrate in Eq. (54) only up to a value Cmax of C, at which the 
photon collision rate becomes equal to the expansion rate, and set 

The exponential factor (with i^jlksTi^ = 25.5) is so sharply peaked at 
C = (max that wc can approximate - - 2Cmax(Cmax - C) in the 
exponent and set equal to Cmax everywhere else in the integral, giving 

Furthermore, Cmax is very close to unity. (For instance, for Um/^b = 7.5, 
we have Cmax = 0.96. That is, we carry the damping integral down to a 
temperature Tl/ Cmax — 3360°K instead of 3100°K.) Hence in this result we 
may as well replace Cmax with unity, so that 

2 1^15(1+0 2(1 +C)V V A J 

This approximation leads to the additional simplification that the damping 
length is independent of most of the parameters appearing in Eq. (56), 
including the ratio Ub/^m- 

There is a smaller additional contribution from the averaging over oscil- 
latory terms, given by Eq. (85). To evaluate this Landau damping term, we 
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will need the ratio At/tL- We noted in Section II that the probability that a 
photon will not be scattered again rises from 2% at about 3400° K to 98% at 
about 2800°K, with very little dependence on any cosmological parameters. 
Matching this to the probabilities calculated from the approximation that 
the probability of scattering in a time interval from t to t + At is a Gaussian 
((it/TT At) exp(-(t - ti)2/At2), and using the relation T oc t~^/^, we find 
At/tL = 0.10, so the contribution to P£ in Eq. (85) has a value 0.0025djj. 
Adding this to the quantity we have calculated from the integral (54) gives 
the total squared damping length 

dl = Vl + AVl 0.029 tl i^^^^lYTl) + 2(ltoO ^ '^'^ ■ ^^^^ 

For instance, for n^/i^ = o.02 (so that ^ = 0.54) Eq. (75) gives du = 
1.61 ti, so = 0.0071 Hence at the first Doppler peak the argument 
of the damping exponential is d'j^k'l ~ 0.07. (This depends very little on ,^.) 
We see that damping is not important at the first Dopper peak, in agreement 
with more accurate computer calculations[17], but is quite significant at the 
second Doppler peak. One effect of damping is to shift the second and higher 
Doppler peaks to lower values of k and i. 

In deriving the wave numbers (86) of the Doppler peaks we also neglected 
the terms proportional to ^/ k^tj^ in the second line of Eq. (83) . At the first 
Doppler peak this quantity is given by Eq. (75) as 



kftl 7r2 

This is 0.20 for Q,Bh^ = 0.03, for which ^ = 0.81, and less for smaller values 
of ^sh^- This approximation is thus fair at the first Doppler peak, and 
becomes excellent at the higher Doppler peaks. 

Finally, we must ask what values of k are small enough so that we can 
ignore acoustic oscillations during the era when the photon energy density 
exceeded the dark matter plus baryon density, during which our analysis 
does not apply. During this era the Robertson- Walker scale factor a{t) went 
as t^^"^, and the speed of sound was l/\/3, so the phase change of acoustic 
oscillations up to the time tc of the crossover from radiation dominance to 
matter dominance was 



Aip = q [ 
Jo 



*c dt 2qtc 2ktL (tca{tL 



\/3a(t) V3a(tc) ^3 \tLa{tc) 
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The redshift zc at the crossover is given by 1+zc = ^m/^-^ = 4x lO^r^M^^- 
During the period from this crossover to the present the scale factor a{t) went 
as t^f^, so the ratio in parentheses is 



tca{tL) ^ j l + ZL _ 1 
tiaitc) ~\l l + zc~ 6.0 VTW? ' 

Using this and Eq. (75) gives 

0.35 / k\ y/l 

For instance, if wc take JIm/i^ = 0.15 and Q^jh^ = 0.03, then Aip ^ 1 at the 
first Doppler peak, indicating that oscillations in the radiation-dominated 
era are becoming important at the first Doppler peak. This is not to say 
that we are making an error of order unity in the argument ip of the sines 
and cosines in Eqs. (83) and (84), but rather that the evolution of the 
perturbations during this much of their oscillations has not been reliably 
calculated. This source of error is mitigated in reference [1] by including the 
effects of photon and neutrino energies on a{t) in calculating the horizon 
distance. 

Our formula (84) for the dipole form factor G{k) raises the possibility 
of a maximum in G{k) at kdu = 7r/2, yielding a "zeroth Doppler peak," 
produced (as the first Doppler peak is not) by the Doppler effect. For 
Ush^ = 0.03 the wave number at this supposed peak is too small for us 
to trust the WKB approximation used to derive Eq. (84) at this peak, but 
the peak in G{k) at kdn = 7r/2 would definitely be there for much smaller 
values of JIb/i^. In particular, the calculations of reference [1] show such a 
zeroth Doppler peak in Q at £ ~ OAbdA/dn for = 0. 

APPENDIX: RELATIVISTIC CALCULATION OF THE 
SACHS-WOLFE AND DOPPLER EFFECTS 



In Section IV we gave a derivation of the Sachs-Wolfe and Doppler ef- 
fects, using heuristic arguments to supplement relativistic results. For com- 
pleteness, this Appendix will present a thouroughly relativistic derivation in 
synchronous gauge, taking into account the possible presence of a vacuum 
energy, which may or may not be constant. This goes over familiar ground, 
first considered by Sachs and Wolfe [18], but I as far as I know there is no 
published treatment of the 'integrated Sachs-Wolfe effect' in synchronous 
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gauge that goes explicitly and analytically into the details presented here, 
including the possibility of a varying vacuum energy. 

A light ray travelling toward the center of the Robertson-Walker coor- 
dinate system from the direction h will have a co-moving radial coordinate 
r related to t by 

= g,,v tot^ydx^'dx" = -dt^ + (a^it) + hrr{rh, t)) dr^ , (91) 

or in other words 

dr / 9 , \ -1/2 1 hrr / X 

The first-order solution is 

1 /■* dt' 

r{t) = s{t) + - J^^ ^hrr is{t')h, t') , (93) 

where s{t) is the zero-th order solution for the radial coordinate which has 
the value tl at i = tL- 

/■* dt' 

s{t) =rL- ^ . (94) 

In particular, if the ray reaches r = at a time to; then 

1 dt 

= s{to) + - / -^Kr is{t)n, t) . (95) 

A time interval 6tL between successive light wave crests at the time 
of last scattering produces a time interval Sto at to given by the variation of 
Eq. (95): 



= dtL 



1 IhrrirLfi-lL) ^ 1 /"*" dt (dh,,\rh,t) 



a{tL) 2 a^{tL) 2a{tL) Jti. a^t) \ dr Jr=s(t)_ 

(96) 



r=rL 



1 ^l/i,,(0,to) 



a{to) 2 a^to) 



(The velocity potential term on the right-hand side arises from the pressure- 
induced change with time of the radial coordinate of the light source in 
Eq. (95).) The total rate of change of the quantity hrr{s{t)n,t)/a'^{t) in 
Eq. (96) is 

d hrr(s(t)h,t) f d hrr{rn,t)\ 1 fdhrr{rh,t) 



dt a^{t) \dt a? it) J ^^.^t) a^{t) \ dr J ^^.^t) 
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so Eq. (96) may be written 

1 1 /irr(0,io) 



= StL 



+ 



a{tL) 2a?{tQ)a{tL) 2a{tL 



to 



dt 



d /hrr{rh,t) 



+StL 



du{rh, t) 
dr 



+ Sto 



r=rL 



1 /irr(0,to) 



a(io) 2 a3(io) 



r=s{t)_ 

(97) 



Hence to first order the ratio of the received and emitted frequencies is 



StL _ a{tL) 
Sto a (to) 



1 f^'^ { d ( hrr{rh,t) 



- a{ti 



/ du{rh, t) 



to 



r=rL 



(98) 



This gives a fractional shift in the radiation temperature observed at time to 
coming from direction n, from its unperturbed value: Tq = TLa{tL) / a{to)- 



AT{h) 



^0 



SW,Dop a{tL)uL/a{to) 



tL 



, , d_f h„-{rh,t) 



r=s{t) 



- 1 



a{ti 



du{rLn, tL) 
dr 



. (99) 



r=rL 



Now we have to think about how to relate the rr component of the 
metric perturbation to the field V' appearing in Section III. In general, the 
metric perturbation may be written as 



h 



The quantity entering into the integrand in Eq. (99) is then 



di V a2(t) 



a(rn, t) + 



d'^f3{rh,t) 

Qj.2 



where 



a 



d_ I A_ 

The field ^ defined by Eq. (23) is given by 

= 3a + 



dt i,2a?) 



(100) 

(101) 
(102) 

(103) 
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We also need a relation between a and /3, which can be taken from the field 
equation for the full metric perturbation [19]: 



■jk , 9 hkk 
+ 



dxWx^ dx'^dx^ dx'^dx^ 
—a^hij + act {hij — S^jkhk^ + 2a^6ijhkk + ^aahij 

= -SttG {Sq - Sp) a^5ij . (104) 

(For simplicity we are here taking the universe to be spatially fiat, which is 
certainly a good approximation at high redshifts, and seems to be a good 
approximation even at present.) The d"^ /dx'''dx^ terms in Eq. (104) give 



A = a B - aaB - 2aaB 



d 



a 



dt V dt 



In terms of the quantities defined by Eq. (102), this is 

a /I d 



a = 



at \aat 



(105) 



(106) 



Hence for a given gravitational potential ip, we can calculate (3 by solving 
Eq. (103): 

i Hi ^^'''-^ 

and then use Eq. (106) to find a. 

Now we return to the fractional temperature shift (99). Using Eqs. (101) 
and (106) lets us write this as 



/AT{n) 
V T 



SW, Dop 



tL 



Qj.2 

1 a 



a{ti 



( au{rLn,tL) 



r=s{t) 



ar 



ItL \9t \adt 
To do the first integral here we note that 

[ a'^Pirh,t)\ 



{a\t) (3irh,t)) 



s{t) 



r=rL 



(108) 



Qj.2 



d_ 
di 



sit) 

a\tf-^^ + a{t)aitmrn^ t) + ait) 



+ a\t) 



a'^(3{rh,t) 

at^ 



at 



+ 3a{t)d{t) 



ar 



r=s{t) 



ap{rh,t) 

at 



+ 



(a{t)d{t) + a^{t)^ (3{rh,t) 



r=s{t) 



(109) 
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The fractional temperature fluctuation (108) may therefore be written 
/ Ar(n) \ _ / Ar(n) \ ^ / AT(n) \ ^ / Ar(n) 

\ T J g-^^ \ T J gai-iy V 2^ / late V 2^ / integrated 

(110) 

where 

fAT{n)\ 2/. ^ (dl3{rLn,t) 



\ ^ _^2(^^) /^5^(^M) -a(i^)d(i^)/3(r^fi,t,.) 

/early V / t=tr. 



early \ / t=t^ 

^ f dp{rh,tL) \ du{rLn,tL) \ , . 



= " (*o) — TT, — + a(to)a(to)/3(0, to) 



V T /late V A=to 



f dfj{rh,to) ^ 

+ "(*o)(^ ) ^ ^ 



r=0 



j^ Ar(n) 



integrated 



-2 ' U^(0 ^^^'^^ +4a(t)d(t) '^^y'^^ +2[a{tyd{t)+a\t)) P{rh,t) 



In evaluating these three contributions to the temperature fluctuation, it 
is helpful to note a relation between (3 and the conventionally defined New- 
tonian potential (p that applies not only for a gravitational field dominated 
by cold dark matter, but also in the presence of a constant vacuum energy. 
Combining Eqs. (26) and (27) gives 

» 2 



r=s{t) 

(113) 



a^ip] =tp . (114) 



dt \4TrGa^pDdt 

Taking into account the relation po on a^"^, an elementary manipulation 
then gives 



d fl d 3 , 



A ^ d f a 



aV (115) 



The equations of the Friedmann model give 

:-47rG(p + p). (116) 



d / d 
dt \a 
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A constant vacuum energy density pv is associated with a pressure pv = 
—pv, while cold dark matter by definition has zero pressure, so as long as the 
gravitational field is dominated by cold dark matter and a constant vacuum 
energy, the right-hand side of Eq. (116) is — 47rGpy , and Eq. (115) then gives 

»^'V^)=o. (117) 



dt \adt 

Comparing Eq. (106) with Eq. (110), we now see that Eqs. (107) has the 
solution 

V^P = i/j (118) 

More specifically, if we define a Newtonian gravitational potential by Pois- 
son's equation 

a-2v^(/) = AttGSqd , (119) 
then Eqs. (26) and (118) show that the Newtonian potential is 

= -|(a2/?). (120) 

This result is not applicable if the gravitational field receives significant 
cotnributions from a varying vacuum energy, but even in quintessence the- 
ories it is reasonable to assume that a vacuum energy density of any sort 
is negligible at and near the time of last scattering. (It certainly must be 
much less than the radiation energy density at the time of cosmological nu- 
cleosynthesis, in order to avoid the production of too much helium.) We 
have also been relying here on the approximation that the radiation energy 
density is much less than the dark matter density at around the time of last 
scattering. Therefore the early-time contribution (111) to the temperature 
fluctuation can be calculated using the relation (118) and tp oc t~^^^, which 
give (3 oc t^^/^. Since here a oc t^/'^, Eq. (120) then gives fS = —t(j)/o? , with 
(j) time-independent. The early-time contribution (111) to the temperature 
fluctuation may therefore be expressed as 

V T /early 3 a{tL) \ Or J r=rL V / r=rL 

(121) 

This yields the Sachs Wolfe temperature shift (68) and the gravitationally- 
induced Dopplcr shift (77) (aside from the terms arising from r = 0, about 
which more later), as well as the pressure- induced Doppler shift (72). The 
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famous factor 1/3 in the first term on the right-hand side arises in 'New- 
tonian gauge' as the sum of a gravitational redshift equal to and a term 
in the intrinsic temperature fluctuation equal to — 2(/)/3, while in the syn- 
chronous gauge used here this term is due entirely to the metric pertur- 
bation. It is a curious feature of synchronous gauge that what we have 
called the gravitationally-induced Doppler shift also arises from the metric 
perturbation. 

It is not appropriate to neglect the vacuum energy at t = so it cannot 
be ignored in the early-time contribution (112) to the temperature fluctu- 
ation. Therefore in general this contribution is not the same as the r = 
terms in Eqs. (68) and (77). Nevertheless, the terms in the early-time con- 
tribution to the temperature fluctuation are only of zeroth and first order 
in h (like the r = terms in Eqs. (68) and (77)) so these terms can only 
affect the multipole coefficients for £ = and £ = 1. 

This leaves the integrated term (113) as the only correction to the results 
of Section IV for i > 2. The integrand vanishes if we ignore the vacuum 
energy and radiation energy, in which case a oc 

i2/3 and /3 oc t''^/^, so the 
integral receives a contribution only for t near to, and is therefore expected 
to be a small correction [20]. Furthermore, although this integral is fairly 
complicated, it has a simple dependence on h. In the presence of a vacuum 
energy, ip(-x,t) can have a fairly complicated dependence on time, but, with- 
out pressure forces acting on the dark matter, its x dependence is the same 
as we found in the absence of vacuum energy, given by Eqs. (31) and (70) 
as 

V'(x,i) = f{t) J d^k e^'^-^k^ek (122) 

with f{t) not proportional to t~^^^ where the vacuum energy is appreciable. 
For a constant vacuum energy, P is then given by Eq. (118) as 

P{^,t) = -f{t) J d^ke^'^-^e^, (123) 

The "integrated" contribution (113) to the temperature fluctuation then 
takes the form 

(^) =2[d^k f'dt^--^^e. 

V / integrated J J to 

X [a\t)f{t) + 4a{t)d{t)f{t) + 2[a{t)d{t) + d'(t))/(t)) • (124) 

It can be shown that this makes an additive contribution to £{i + l)Ci that 
for large £ goes as 1/^, with no interference between this contribution to the 
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temperature fluctuation and the other contributions. [21] For a time- varying 
(but spatially constant) vacuum energy the function /3(x, i) does not satisfy 
the relations (118) and (123), but Eq. (107) shows that its spatial Fourier 
transform is nevertheless just proportional to ek for large k, so the integrated 
term still makes a contribution to i{i + 1)C^ that is proportional to 1/^ for 
large i. 
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